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The Bose-Einstein condensation in the hard-core boson limit (HCB) of the Bose-Hubbard model with two local 
^ ', states and the particle hopping in the excited band only is investigated. For the purpose of consideringthe non- 

ergodicity, a single-particle spectral density is calculated in the random phase approximation by means of the 
bJQj temperature boson Green functions. The non-ergodic contribution to the momentum distribution function of 

particles (connected with the static density fluctuations) increases significantly and becomes comparable with 
the ergodic contribution in the superfluid phase near the tricritical point. 
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1. Introduction 

The Bose-Hubbard model (BHM) yj] is widely used in describing the thermodynamics and the dy- 
namics of Bose atoms in optical lattices |2, 3]. Such systems demonstrate the phase transition into the 
phase with Bose-Einstein (BE) condensate at very low temperatures. Thus, the system can be in the nor- 
mal (NO) phase (the state of the so-called Mott insulator, MI), or in the phase with the BE condensate (the 
superfluid state, SF). The BHM can be used for description of various phenomena such as the quantum 
derealization of hydrogen atoms adsorbed on the surface of transition metals |4], quantum diffusion of 
light particles on the surface or in the bulk |5], thermodynamics of the impurity ion intercalation into 
semiconductors @. 

The BHM has been intensively studied for the two recent decades starting with the pioneering work 
0], where the thermodynamics of the model was considered in the mean field approximation (MFA) for 
a particle hopping with an exact treatment of the Hubbard repulsion. Such a description revealed the 
first order NO-SF phase transition for particles residing in the ground state of local wells. A prospective 
generalization of the BHM consists in taking into account the excited vibrational states of bosons in local 
potential minima of lattice sites. For example, intersite particle hopping through the excited states is 
much easier in the case of quantum derealization or diffusion |8, 9] as well as in the case of optical 
lattices 1 10]. However, a possible BE condensation in the excited states was considered for the case of 
optical pumping in order to increase their occupations 1 11]. In this case, the orbital degeneration of the 
excited p-state can lead to the appearance of a special type of the BE condensate. 

A thorough study of the phase transition into the phase with BE condensate in the BHM with two 
local states and the boson hopping only in the excited band was carried out in works jl^El]- In the MFA 
and the hard-core boson (HCB) limit, the instability related to the NO-SF transition was studied which 
occurs at excitation energies lower than the absolute value of the hopping parameter. Conditions of the 
change of the above mentioned phase transition from the second order to the first one were studied. The 
respective phase diagrams (0,p) and (I^L/i) were analyzed. A possible phase separation into NO and SF 
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phases at a fixed concentration was revealed. Two-time boson Green's functions (GF) and a single-particle 
spectral density were calculated in the random phase approximation (RPA). The structure of excitation 
spectra of the particle and hole types in the region of the NO-SF phase transition was investigated. 

Some physical many-particle systems do not hold the ergodicity hypotesis. It means that there are 
regions of the phase space of a system which cannot be reached by the trajectory of the system evo- 
lution. From practical viewpoint, non-ergodicity leads to a distinction between isothermal and isolated 
responses (susceptibilities) of a system caused by the zero-frequency term present in the isothermal re- 
sponse only. In this case, the use of the "isolated" (Kubo) response leads to erroneous values of such 
physical characteristics as compressibility, specific heat and susceptibility for the systems described by 
means of the Ising model (3, EH (a zero Kubo response), the pseudospin-electron model (3 an( i the 
Falicov-Kimball model (the non-ergodic contribution describes the singularity related to the phase tran- 
sition (lJEll)- The BHM in the HCB limit, as will be shown below, also belongs to the non-ergodic class. 
The mentioned problems can be solved using the temperature (Matsubara's) GF 1 19] (see, e.g., a scheme 
for calculation of many-time correlation functions 1 20]). 

In the present article we study the non-ergodic contribution to the particle momentum distribution 
function by the example of the two-state Bose-Hubbard model in the HCB limit. A detailed investigation 
is performed in the region close to the tricritical point where the order of the phase transition to the SF 
phase changes to the first one. 



2. Thermodynamics of the model in the hard-core boson limit 

The Hamiltonian of the BHM takes into account the tunnel hopping of particles between the nearest 
neighbour lattice sites and the mutual repulsion of the particles located in the same potential well fl]: 

& = E r '7 b i h i + 7 E n i («i - 1) - ME rH ■ (2.1) 

ij Z i i 

Here, t%j characterizes the hopping energy, U is the energy of the Hubbard intrasite interaction (U > 0), 
\i is the chemical potential of the boson particles, bt,bj are the Bose operators of particle creation and 
annihilation on the z-th site in the ground vibrational state. Taking into account the first excited state and 
considering the particle hopping only between these excited states, one can generalize the Hamiltonian 
(ID as follows: 

A = (e - J>+ ^ + (£' - M) E c t * + ^ E A (n? - 1) 

i i i 

+ ^7 !>?(»/- 1) + tf* !>?»?+ E t'ij c+ i c i - ( 2 - 2 ) 

i i ij 

where ct, a are the Bose operators of the particles in the excited state, e {e') is the energy of the particle 
in the ground (excited) state, Ut, U c , Ubc are the parameters of the Hubbard repulsion. 

Starting from the basis of the single-site states \ i; m, np, formed by the particle occupation numbers 
(eigenvalues of the operators rfi - btbi and = etc;), one can introduce the Hubbard operators jilll 

X^ m ' n '' m ' = \i;n,m){i;n' ,m'\. (2.3) 

Let us restrict our study to the HCB case (when n + m =£ 1). In the considered model, such a situation 
is reached at Ub,U c ,Ub C — - oo. In this limit, the single-site Hamiltonian becomes three-level with the 
respective energies Xq — 0, X\ — X2-8- p. (the following shortened designations for single-site states 
are used: |0) = |00>, |1) = |10), \2) = |01». Here 8 = e' - e is the energy of transition to the local excited 
vibrational state. Respectively, 

b i = X m - w = Xf\ Ci = X 00fil =X° 2 ; 
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Finally, the Hamiltonian {2.2) is formulated in the X -operator representation 



(2.4) 



The BE condensation takes place in the band that originated from the hopping between the excited 
states of the neighbouring wells (see below). Thus, the order parameter is equal to the average of boson 
creation or annihilation operators £ = (Xf°) = (X? 2 ) = (ct) = (c,-». 

Let us separate a MFA part of the Hamiltonian H2.4) 



(2.5) 



ip i 



where t' Q is the Fourier transform of the hopping energy at q — (hereinafter we consider t' Q < 0, 
i.e., assuming symmetrical excited states, while t' Q > for antisymmetrical ones |22j). The £ parameter is 
obtained from the self-consistency equation £=Z _1 Sp[X 20 exp(-pHMp)], where Z = Spexp(-/3/ZMF)- 

This approximation corresponds to the exact solution for the case of an infinite-range hopping. A 
strict derivation of the above statement for the standard BHM has been given in the work 12311 (see also 
12411 ) based on the approach in the style of the Bogolyubov-Tyablikov variational method. 

Finally, a complete Hamiltonian looks as follows: 



H = H M f + E A) ( x f° - o ( X T - o. 

The MF Hamiltonian Hmf can be reduced to a diagonal form by the transformation 



(2.6) 



( l0> 1 




' cosd 





-sin£ 


ID 







1 





I 12) J 




k sin d 





COS d 




where 



cos 20 = (<5 -/*)/£, sin2# = 2|f K/£'; E= yj{8-y,) 2 +4(t' O 2 - 
In terms of operators X^ v = |jtt)<v| 



(2.7) 



(2.8) 



(2.9) 



Energies of single-site states 1^ in the phase with a broken symmetry (£ + 0) are equal to Ao,2 = \ {S-fi+E), 
Xi = -ji. 

It is convenient to introduce the following linear combinations 



X 22 ) 



■\a y .=X\ 



02 



20 



(2.10) 



Operators a" have the same properties as the spin operators for the case S = 1/2 and fulfil the same 
commutation rules (with an important difference: the anticommutator of ert and cr7 operators is equal 
to X?° + Xf 2 instead of unity as in the ordinary case). 

Thus, the X-operators are represented on a new basis as follows: 



X 20 = sin2tfof + cos2#crJ 
with the MF Hamiltonian in the form 



-ler. 



8 + . 



X 



02 . 



(2.11) 



(2.12) 
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The Hamiltonian i'Z.G) expressed in terms of 
operators <2.10> is a three-state generalization of 
the standard Hamiltonian of the HCB model 12511 . 
The presence of the third level |1> does not effect 
the pseudospin dynamics but modifies the occu- 
pations of levels |0) and |2), changing this way the 
thermodynamic behaviour of the model. More- 
over, the application of pseudospin formalism al- 
lows one to utilize the usual scheme of the random 
phase approximation (RPA) for Green's functions 
of spin models, originated from the well-known 
Tyablikov decoupling for the Heisenberg model 

El. 

Let us make a short revision of the thermody- 
namics of a three-state model in the MFA. Using 
the Gibbs distribution with the Hamiltonian <2.9> 
and taking into account that the average values of 
transverse pseudospin projections (o4) are equal 
to zero, one can obtain the equation for the order 



parameter f = (X? 2 ): 
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Figure 1. Phase diagrams (Q,/j) at various values of 
S (I £g I = 1). The solid line corresponds to the second 
order phase transition and the dashed line denotes 
the first order transition. 



f = Z~ l E 



-PA.0 _ 



)■ 



(2.13) 



The case of £ = corresponds to the NO phase. One can obtain nonzero solutions (indicating the appear- 
ance of the BE condensate) from the following equation 



Z~ l E 



V l(e 



-/M : 



(2.14) 



It has been shown in work 1 12] that in the region \i < 0, nonzero values of £ gradually appear at 
the second order transition while in the region [i > at a low enough temperature, the curve £(/i) is 
S-shaped which indicates the first order phase transition. The line of this transition is defined by the 
equality of values of the grand canonical potentials Q, = -QlnZ in phases NO and SF. The change of the 
phase transition order and the localization of the respective tricritical points are presented in the phase 
diagrams in figure [T] 

Unlike the standard (two-level) HCB model (where the Bose-particles remain in the local ground state, 
the excited local states are neglected and phase transitions between NO and SF phases are always of the 
second order), the phase diagrams in figure[l]are asymmetrical. 



3. Temperature Green's functions and the non-ergodic contribution 

The spectrum of boson particles and the single-site spectral density for this model are calculated in 
the RPA using the two-time GF in work 1 13] . Herein below we will perform similar calculations in the RPA 
but applying the temperature GF. 

Assuming the hopping between the excited states only, the dynamics of the model corresponds to the 
HCB model and we consider only GF constructed on the operators which describe the transition between 
the states |0> and \2) 

( T T X 20 {j)Xf (t') ) = f + ( T r (Xf - S) T {X] 2 -Or'), 

(T T X 20 X° 2 ) q ^ n = pe6(q)6(co n ) + (TAX 20 -0(X° 2 -V) qb)n . (3.1) 
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Using the above mentioned "pseudospin" representation H2.10K one can express (3.1) in the form 



( T r {X 20 - {X 02 -0) at0 =^ zz sin 2 2-0 + -<^ ++ (cos 2 2d - 1) 



^ +_ (cos2# + l) 2 + |^- + (cos2fi- l) 2 + i^"-(cos z 25- 1) 
l 



.„, ; + ^ z -)(cos2fl+l)sin25 
-i(^" z + ^ z+ Hcos2#-l)sin2# ( (3.2) 

a " = ( T r {a a - {<J a )) T {<T p ~ (o p )) T ) . (3.3) 
Temperature GF in the RPA are derived by means of the Larkin equation 

VP = ^ + i 1 * 8 ^ , (3.4) 

yd 

where unperturbed GF are expressed through the single-site GF and the average values of the pseudospin 
operator 

% + ~ = -2g°(co n )((T z )o, % z = pb'8{oj n ), 

% + = -2g\-b) n ){<j z ) , sV») = k^I> ( 3 - 5 ) 

KPE) = {<j z ) Q = | tanh ^ , b ' = Wm= < CT '>o ~ < CT *>o > 

and the parameters of the effective interaction look as follows: 

7 ++ (?) = j (<7) = \ (cos 2 2d -\)t' q , 
j + ~ iq) = j " + (?) = | (cos 2 2d + 1) ^ , 
j zz {q) = 2t' q sm 2 2d, 

j +z {q) = j z+ (q) = f z (q) = j z ~{q) = -^sin20cos20. (3.6) 

The equation of the type <3.4) for the ordinary HCB model was derived in work l27ll . 

The momentum distribution of particles in the excited state is obtained as the sum over Matsubara's 
frequencies 

n q - £ 2 8{q) = - - £ e 1 *"^ ( T T {X 20 - {X 02 -Q) ^ (3.7) 

P U)„ 

from the temperature GF 

'cos2^-O fl cos2# + 



( TV (X 20 - {X 02 - $)) q(i)n = (ct z ) | — ^ + - ; - q \+<3i{q)-6{<o n ), (3.8) 



<b q = — [£(cos 2 25 + 1) +4<(7 z >^cos 2 2#] , (3.9) 

£,7= [(£ + 2<(T z >^](£' + 2(cr z >^cos 2 2^]] 1/2 . (3.10) 

An explicit form of the spectrum depends on the phase where the system resides. For example, for NO 
phase 

£ (NO) _ g _ ^ + 2 ( a z ) t' q , E = 8-n, (3.11) 
while there are two branches +£„ F) for the SF phase 



e r= 2 i<^>i{(i^i+4)^oi+44pw 



1/2 



E = 2\t' \{a z ). (3.12) 
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Figure 2. Dependences (<j z )q and b' on chemical potential at different temperatures (5 = 0.1, | t' Q | = 1). 



This result formally coincides with its analogue for ordinary (two-state) HCB model (which can be reached 
in the limit S — ► -00). However, an important difference lies in the fact that the values of (er z ) and d are 
defined by an equation for the three-state system 12.12) . 
There are three parts in the momentum distribution: 



77„ = eS(q) + <o - : >„ |0,coth^ -cos2tf) + jWq) ■ 



(3.13) 



BE-condensate 



single-particle 
excitations 



The first and the second terms correspond to the results obtained in work 1 13] while the non-ergodic part 
(related to the static density fluctuations) can be obtained only by means of temperature GF: 



2/sin 2 2tf 



(l+2^ t ' q cos 2 2#) (l+2^% cos 2 2tf +2^' t' q sin 2 2#) ' 
1 ~ 



{X 00 + X 22 ) - {a z ) 



Z\2 



(3.14) 



The factor sin 2d is non-zero in the SF phase only, so the existence of the non-ergodic part is also limited 
to this phase: 

| = 0, NO phase, 



+ 0, SF phase. 



Due to a sophisticated structure of the nonergodic part, a thorough analysis of its behavior is rather 
complicated. Dependencies of {a z ) and b' parameters, entering the formula H3.14) , on chemical potential 
and temperature are presented in figure [2] As one can conclude, the non-ergodic contribution is "frozen 
out" at low temperatures (because the hopping is possible between the excited states only). 

The momentum distribution, including according 13.13) all three parts, rapidly increases approaching 
the centre of the Brillouin zone (BE-condensation). The relative role of the non-ergodic part in the total 
momentum distribution is illustrated in figures[3H5] In the SF phase E-2\ t' a \{a z )a, therefore 



l + 2 ^£ f 'cos 2 2^ 
E * 

1 



q^0 



-<Si{q)^oo 

near the line of the second order SF-NO transition. Enhancement of the non-ergodic contribution in the 
vicinity of the tricritical point is connected with the increase of the particle number fluctuations (the 
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latter are characterized by the derivative dn/d/i, that reduces to the b' parameter in the simplest approx- 
imation). It is confirmed by the behaviour of the dependence at different temperatures obtained and 
illustrated in flil . 

The total occupation of the excited state can be calculated by the summing over the wave vector (here 
we substitute the summation by the integration with the model density of states go (z) for a simple cubic 
lattice Q) 

"rpa- ^£"4 = Jn z g iz)dz. (3.15) 

It is interesting to compare partial contributions of the BE-condensate, the single-particle excitations 
and the non-ergodic part in the total occupation of the excited state: 

"rpa = » - ? + "sp + "ne • 

The most pronounced non-ergodic contribution is visible near the line of the second order phase transi- 
tion in the vicinity of the tricritical point (figure[6}. 




Chemical potential ^ Chemical potential ^ 

Figure 6. Comparison of the contributions to the excited state: the BE-condensate, the single-site contri- 
bution and the non-ergodic part (8 = 0.1, | tL | = 1). 



4. Conclusions 

The non-ergodic part of the momentum distribution (related to the static density fluctuations) in the 
excited state for the BHM in the HCB limit is calculated by means of the temperature GF. This part be- 
comes nonzero only in the SF phase. Its value rapidly grows and becomes on par with the ergodic part in 
the SF phase near the tricritical point (where dnld\i — * 00). 

The questions remain, is the non-ergodicity (and the presence of corresponding contributions to the 
particle momentum distribution) a characteristic feature of the hard-core boson limit only and what 
could be at the transition to the standard BHM. Such a transition takes place when a constraint of occupa- 
tion number n,- is removed and the on-site repulsion energy U becomes finite. When the non-ergodicity 
would disappear in this case, it could mean that the effect is caused by the change of the particle statistics 
(due to the transition from Pauli- to Bose-operators). This problem needs a separate and more detailed 
consideration. 

Dynamical characteristics of the model obtained in the present study and in the work (l3J] are gen- 
erally based on the average values calculated in the MFA. There is a slight discrepancy between these 
characteristics and the respective results in RPA. For a complete self-consistency, one should use the ap- 
propriate expressions for averages, starting with the grand canonical potential obtained by summation 
of diagram series consisting closed cycles. 
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AsocTaHOBa MOfle/ib Eo3e-Xa66apfla b rpaHi/uji >kopctki/ix 
6o30hIb: HeeproflMHHicTb Ta 6o3e-i<OHfleHcai^in 

LB. CrachOK, O.B. Be/ii/iHKO 

IhtCTHTyT (J>i3HKl<l KOHfleHCOBaHHX CWCTeM HAH YKpaiHH, 79011 /lbBiB, By/1. CBCHL4iL4bKOrO, 1 

flooiiflxyeTbCfl 6o3e-KOHfleHcai4m b rpaHWL(i xopcTKWx 6o30HiB y MO,qe/ii Bo3e-Xa66apfla 3 flBOMa /lOKa/ibHUMn 
CTaHaMU npw nepeHOci 6o30HiB /lnwe y 36yflxeiHM 30Hi. 3 MeTOio BpaxyBam-ifl HeeproflWHHOCTi 0flH04acTWH- 
KOBy cneicrpa/ibHy rycTHHy OTpuMaHO b Ha6/in>xeHHi xaoTH4Hwx (|>a3 3a flonoMoroio TeMnepaTypHHX 6o30hhwx 
cpyHKi^ii/i TpiHa. HeeproflH4HwR BHecoKflO cpyH KL^iT po3nofli/iy nacTHHOK 3a iMny^bcoM (noB'fl3aHwti 3i CTaTUHHM- 
mh cfi/iyKTyaLjiflMW rycTHHn) cyTTCBO Hapodae i crae cniBMipHUM 3 eproflM4HOK) nacTHHOK) b Hafln^nHHiR cba3i 
6\nn TpuKpnTUHHOi tohkh. 

K/iioMOBi c/ioBa: MOfle/ib Bo3e-Xa66apfla, >KopcTKi 6o3ohi/i, 6o3e-efit-mjTeiiHiBCbKa KOHfl,eHcau,iR, 36yfl>KeHa 
30Ha, Heeprofln^HicTb 
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